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We investigate the classical evolution of a <f> 4 scalar field theory, using in the initial state random 
field configurations possessing a fractal measure expressed by a non-integer mass dimension. These 
configurations resemble the equilibrium state of a critical scalar condensate. The measures of the 
initial fractal behavior vary in time following the mean field motion. We show that the remnants of 
the original fractal geometry survive and leave an imprint in the system time averaged observables, 
even for large times compared to the approximate oscillation period of the mean field, determined by 
the model parameters. This behavior becomes more transparent in the evolution of a deterministic 
Cantor-like scalar field configuration. We extend our study to the case of two interacting scalar 
| fields, and we find qualitatively similar results. Therefore, our analysis indicates that the geometrical 

\^ , properties of a critical system initially at equilibrium could sustain for several periods of the field 

oscillations in the phase of non-equilibrium evolution. 
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I. INTRODUCTION 
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The classical dynamics of scalar field theory has been extensively studied in the literature. In most investigations 
a lattice discretization of the scalar field is used, reducing the problem to the study of the dynamics of a system of 
non-linear coupled oscillators [l[ . The central question in this case is the evolution of the system towards a thermalized 
stationary state. In the early days Fermi, Pasta and Ulam [2j have obtained deviations, even for large times, from 
the naively expected equipartition of the energy among the different oscillators. Through the efforts to explain these 
results, it became clear that, for appropriate initial conditions, a variety of periodic solutions (breathers, solitary 
43 ' waves) [H, defined on the non-linear lattice, exists. Therefore, the choice of the ensemble of the initial configurations 
strongly influences the long time behavior of the system dynamics. Recent works [H, @ show that for a random 
ensemble of initial configurations a sufficiently large system relaxes to the usual equilibrium distribution, but the 
£S| . corresponding relaxation time strongly depends on the parameters of the theory. 

^ ■ In the present work we reconsider the classical dynamics of a scalar field in 1 + 1 dimensions adopting a different 
point of view: As initial conditions we use an ensemble of scalar field configurations possessing a non-conventional 
profile, inspired by the order parameter fluctuations of a critical system at thermal equilibrium. These configurations 
generate a fractal measure on the lattice characterized by a corresponding fractal mass dimension. We do not consider 
here the dynamical process responsible for the formation of such a critical state, but we concentrate on its evolution 
iy-) ■ once it has been formed. Our aim is to investigate the deformation of the initial fractal measure as the system evolves 
according to the classical equations of motion. In particular, we are interested in determining the time scale for 
which signatures of the initial fractality survive, and leave their imprint in appropriate observables. We find that the 
initial geometry is successively deformed and restored again, with a frequency determined by the field oscillations. 
This behavior seems to be generic, since it is observed for both random as well as deterministic (Cantor-like) fractal 
measure. Moreover, we study in detail the influence of an additional thermalized non-critical scalar field, described 
43 i initially by configurations with conventional geometry, coupled to the system, as well as the dependence of the 
C^' corresponding characteristic time scales on the parameters of the theory. Finally, we discuss the applicability of our 
model to the out-of-equilibrium evolution of an isoscalar condensate formed near the Quantum Chromodynamics 
(QCD) critical point during a heavy-ion collision experiment. 

The paper is organized as follows: Section II is divided in two subsections. In subsection IIA we present our 
model considering a single self- interacting scalar field. In IIB we describe the generation of the ensemble of initial 
configurations for a scalar field corresponding to a random fractal measure on the lattice with a given fractal mass 
dimension. Section III contains the numerical results of the single field case and section IV the corresponding results 
for the evolution of field configurations initially generating a deterministic fractal measure with Cantor-like structure. 
In section V we analyze the dynamics of two coupled scalar fields. Finally, in section VI we summarize our results 
and discuss their relevance to the phenomenology of out-of-equilibrium critical systems. 
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II. THE MODEL 



A. Equations of motion 

Our model consists of a classical scalar field ct obeying the usual 4 -dynamics described by the Lagrangian density: 

£ = \d^a - V(a) (I) 

with the potential 

V(a) = ±(a 2 -l) 2 -Aa, (2) 

where A and A are the coupling parameters of our model. All the quantities (ct, A, A, as well as the space-time variables) 
appearing above are chosen dimensionless. Following the cr-model we assume that the Z<i symmetry (ct — > — ct) is 
broken only through a linear term in the potential, setting the coefficient of the cubic term to zero. Furthermore, we 
have absorbed one more parameter by rescaling the field as well as the space-time units. Thus, only two parameters 
remain in the potential term. We consider the dynamics of the scalar field in f + 1 dimensions. The corresponding 
equation of motion is 

a - ct" + Act 3 - Act - A = (3) 

where dot represents time derivative and prime the spatial one. To proceed numerically we have to discretize eq.© 
on a lattice. This reduces the system to a chain of non-linear coupled oscillators. We use the following leap-frog 
discretization scheme: 

dt 2 

°i = 2(7 i - °i + K+l + a i-l ~ 2(T t ) ~ dt l^ a i > ~ Xa t ~ A \ ' ( 4 ) 

where dx is the lattice spacing, dt is the time step, the upper indices correspond to time steps and the lower indices 
to lattice sites. As usual we perform an initial fourth order Runge-Kutta step to make our algorithm self-starting. 
We are interested in studying the evolution of the above system determined by eq.((3]), using an ensemble of initial 
field configurations possessing a non-conventional profile, characterized by a fractal mass dimension. The motivation 
of this choice and the details of constructing such an ensemble, defined on a f-dimcnsional lattice, are given in the 
next subsection. 



B. Generation of initial ensemble of a-configurations corresponding to a random fractal measure 

The absolute value of the u-field introduced in the previous subsection is interpreted as local density, and the 
corresponding fractal behavior is described by a fractal measure demonstrated in the dependence of the mean " mass" 
on the distance R around a point xq defined by: 

m(x ,R) = ( \a(x - x )\d D x), (5) 
Jr 

obeying the power law 

m(x ,R)~R Df (6) 

for every xq- Df is the fractal mass dimension of the system 0, 0, H| and the mean value is taken with respect to the 
ensemble of the initial cr-configurations. The production of a cr-ensemble possessing the fractal measure described in 
cqs. (I5I6[) . has been accomplished in [9(. It is based on the observation that a scale invariant free energy of the form: 

F[a}= ^d D x{l(Vaf+ga s + 1 } 1 (7) 

when introduced as a weight in the partition function: 



Z= js[a]e-^, 
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of the scalar field a, generates piecewise constant configurations leading to an ensemble possessing fractal mass 
dimension (according to the definition above): 

In the following we will use 6 — 5 and the dimension D = 1, therefore the corresponding fractal mass dimension is 
Df = 5/6. It must be noted that the free energy ([7]) for D = 3, 6 — 5 and g — 2 describes the effective action of the 
3D Ising model at its critical point (Toj . 



In practice, to produce critical configurations on a lattice we use the following algorithm 11|: We perform a 
random partitioning of the lattice in elementary clusters of different size £. Thus, each cluster consists of several 
lattice points. Within each cluster the value of a is assumed to be constant. To obtain the values of a in the different 
clusters we use a uniform random distribution. Each cr-field configuration is weighted by a factor e _r [^_, where T[cr] is 
calculated for the given configuration, using the coupling value g = 2. We use a Metropolis algorithm [12| to perform a 
random walk in the configuration space of the cr-field. We perform firstly « 25000 initial algorithmic steps in order to 
achieve equilibrium. The a ensemble is then formed by recording a large number (~ 10 4 ) of statistically independent 
cr-configurations . 

With this procedure we acquire an ensemble of field configurations generating a random fractal measure on the 
lattice as a statistical property after ensemble averaging. This property is not reflected in the geometry of either 
one configuration neither in their average (o~(x)}, which is depicted in fig.[TJ but is produced only through the entire 
ensemble. Each one of the configurations as well as their average, have a continuous power spectrum (P(f) ~ / — ly/2 ) 
resembling a colored noise profile, and the mean value of the field (spatial average) is almost zero. 




FIG. 1: The a-field on the lattice averaged over the ensemble of the initial configurations. 



The produced ensemble possesses the property ([6]), where now the fractal mass dimension is determined by the 
power-law behavior of m(xo,^) = (J^° + |cr(x)|dx) around a random xq, averaged inside clusters of size £ 1 . The 
(fx° + ^ \o~(x)\dx) versus £ figure is drawn as follows: For a given xo of a specific configuration we find the size £ of the 

cluster in which it belongs and we calculate the integral J*° + ^ \cr(x)\dx, thus acquiring one point in the (f^° \cr(x)\dx) 
vs £ figure. For the same Xq we repeat this procedure until we cover the whole ensemble, and the aforementioned 
figure is formed. Averaging in xq obviously does not alter the results, since m(xo,£,) ~ m(xo + with I spanning 

the entire lattice. In figOJi we observe that in the log-log plot of (f*° + ^ \cr(x)\dx) vs £, the slope ip, i.e the fractal 
mass dimension Df according to ©, is equal to 5/6, which is the theoretical value calculated from |J5]). within an 
error of less than 0.3%. 



1 Note that this fractal mass dimension must not be confused with the fractal dimension of the corresponding curve, which in this case 
is greater than 1. 
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FIG. 2: Potential V{o~) given by flp, for different values of A and A. 
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FIG. 3: Time evolution of the ensemble average of the spatial mean (a) for different values of A and X in dimensionless units. 
Mind the different vertical axis scale in the last graph. 



III. NUMERICAL RESULTS. SINGLE FIELD 

We study the evolution of the system determined by equation ([3]) which we solve in 1-D 2000-site lattice, using 
as initial conditions an ensemble of 10 4 independent cr-configurations on the lattice generated as described above, i.e 
possessing fractal characteristics. The initial time derivatives of the field, that is the kinetic energy, are assumed to 
be zero, since this is a strong requirement of the initial equilibrium. The ensemble population is by far satisfactory 
since the results are independent of it as long as it is larger than 6 x 10 3 (numerically tested), and furthermore they 
are independent from the number of the lattice sites provided that it is greater than ~ 10 3 . We investigate the 
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evolution of m(xo,£){ = (J^° \cr(x)\dx}) which initially is a power law - ^(°), where V(0) = D s = 5/6 (fractal 
mass dimension). 

In order to understand the dynamics of the cr-field within this ensemble we firstly consider the evolution of the field's 
mean value (cr), as well as its standard deviation ((5er) 2 = (a 2 ) — (c) 2 - The averages are taken over all statistically 
independent configuration. We classify the dynamics in six different cases, produced through a suitable choice of the 
parameters A and A, according to the corresponding form of the potential plotted in fig [2] In fig. [3] we depict the 
evolution of (a) for these cases. For zero A and A, the potential term vanishes and the mean value of the field remains 
constant and equal to a very small value determined by the initial conditions. For non-zero A the field oscillates 
around the potential minimum and the oscillation amplitude as long as the frequency increase with A for fixed A. 
This is due to the fact that the minimum value of the potential decreases and at the same time the value of cr, for which 
the potential minimum occurs, increases (see figs. [2] and [3]). Note furthermore that due to the quadric anharmonic 
term in the potential the time mean value of the oscillations is slightly smaller than the potential minimum. On the 
other hand, for small A, (a) oscillates almost harmonically, and the amplitude together with the minimum decrease 
with increasing A, while the frequency increases. For A > 1 the oscillations damp relatively fast (the larger the A 
values the faster the rate, as the potential becomes steeper), leading to a stabilization of (a). Finally, in fig. [4] we 
show the evolution of the ensemble average of the standard deviation of the cr-field for the same parameters as above. 
Note that the standard deviation can be quite large even if (ex) remains small, such as in the A = 1, A = 100 case. 



A=0, X=0 M ] A=1,).=1 2,s ] A=10, i=1 




FIG. 4: The standard deviation time evolution of the ensemble average of the a-field for different values of A and A in 
dimensionless units. 



As a next step we calculate the evolution of m{x ,£) defined previously. As time passes the power-law form of 
m(xo,£) (~ remains, but the corresponding exponent ip(t) increases, approaching the value ip(t) » 1, when the 

signs of the initial fractal geometry disappear and a conventional pattern establishes, as is depicted in fig.[5j However, 
a more detailed analysis for greater time intervals reveals a remarkable phenomenon. In the solid line plots of fig. [5] 
we show the evolution of ip (each ip value coming from a linear fit). It can be clearly observed that the characteristic 
exponent i[i(t) after reaching the value 1, fluctuates and for particular times becomes almost equal to ip{0) = 5/6. 
Thus, the initial fractal measure is restored repeatedly within the considered ensemble. A simultaneous study of this 
graph and fig. [3] is illuminating, since we observe that the reappearance of the fractal mass dimension, describing the 
initial measure, occurs when the mean field value becomes almost zero. 

If we interpret the cr-field as a system of coupled anharmonic oscillators with initial field values around zero and 
zero kinetic energy, we expect an almost simultaneous pass from their turning points, due to synchronization. The 
energy transfer between the different oscillators takes place through the spatial derivative, as can be seen in equation 
([3]). Therefore, if it is small compared to the other terms (which is the case in general), the oscillators do not mix 
significantly, and due to their initial zero kinetic energy they roll to their common minimum, they oscillate around it 
and then return to their starting position almost in phase and with kinetic energy close to the initial one, i.e close to 
zero. So, every time the kinetic energy and the mean field value become zero, the system reaches to a state similar 
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FIG. 5: {f*° +i \o(x)\dx) versus £, for successive times t = 0,t = 0.25, t = 0.5, t = 0.75, t = 1, t = 1.25, /or ,4 = 1 and A = 1 
in dimensionless units. We observe that although the initial fractal mass dimension changes, the power law property remains 
valid as time evolves. 
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FIG. 6: Time evolution of the slope ip (solid line) and its time average {ip)t (dotted line) for varied A and X in dimensionless 
units. 



to the initial one, leading to the restoration of the initial fractal measure. On the other hand, when the oscillators 
roll away from zero, the geometrical characteristics of the ensemble are changed by the dynamics, since they acquire 
large field values and large kinetic energy (compared to the initial ones). 

In figs [6)3, c,d we can clearly observe this behavior. Moreover, ijj(t) in figsEJj and f is also easily interpreted since 
in these cases, following closely the evolution of the a mean value, it stabilizes after some reciprocations. Finally, the 
behavior of tp for A = 0, A = depicted in fig. [6^ is also expected since in this case ([3]) reduces to the simple wave 
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equation, a values stay always around zero, i.e around their initial values, and ip oscillates around a value between 
5/6 and 1, namely the signs of the initial fractal geometry are always visible. The dotted line plots of fig. [5] depict 
the time average of ip defined by {yb) t = \ / ^(f)df. As it can be observed, the influence of the initial fractal mass 
dimension is still visible in this time integrated measure. 

Lastly, in order to have a clearer apprehension of the aforementioned phenomenon, we perform some tests. In fig. 
Wfr and [TJd we present the evolution of the system with random initial conditions, prepared choosing the a value at 
each site from a uniform distribution, for A = 1 and A = 1. The mean field value oscillates around the minimum as 
before, and the slope ip which initially is obviously one, corresponding to a non-fractal system, remains equal to 1 as 
expected, independently of the field motion. Following the arguments referred above, one could say that every time 
(a) returns to zero the system enters in a state similar to the initial one, and vb(t) remains always 1 since the initial 
state is characterized by tp(0) = 1. In figs[71: and [31 we evolve our system using initial conditions corresponding to a 
fractal measure with mass dimension ?/>(0) = 5/6, but with random non-zero (actually quite large) kinetic energy, for 
A = 1, A = 1. In this case the finite value of the initial kinetic energy (different for every oscillator) forbids the return 
of the system to a state close to the initial one, each time the mean value of the field passes through zero, suppressing 
the approach to the initial fractal state. (Note that due to the large initial kinetic energy the system oscillates around 
both minima). 

Additionally, in figs[7]3 and[7f we perform the following scenario: We evolve the initial ensemble, possessing fractal 
mass dimension ip(0) = 5/6, for A = 1, A = 1, introducing by hand a three orders of magnitude larger coefficient to 
the spatial derivative term in the equation of motion. As expected, the enhanced diffusion in this case induces very 
strong mixing of the oscillators, and as a consequence the system never reacquires the initial fractal characteristics 2 . 




t t t 



FIG. 7: Mean field, ip{i) (solid line) and (ip) t (dotted line) evolution for various scenarios, a) and b) correspond to random 
non-fractal initial conditions for A = 1 and A = 1, c) and d) correspond to initial conditions with fractal mass dimension 
t/>(0) = 5/6 and non-zero, randomly chosen, initial kinetic energy, for A = 1 and A = 1, and e) and f) correspond to initial 
conditions with fractal mass dimension ip(0) = 5/6, zero kinetic energy, but three orders of magnitude larger spatial derivative 
term, for A = 1 and A = 1. 



2 Note that the increased diffusion can also be achieved using the modified form of the potential: V(a) = 4 (a 2 — v 2 ) 2 — Aa, taking small 
A and A and large v. 
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IV. EVOLUTION OF A SCALAR FIELD CONFIGURATION WITH DETERMINISTIC FRACTAL 

MEASURE 

In order to acquire a better comprehension of the mechanism of the aforementioned phenomenon, we investigate 
the evolution of one Cantor-like scalar field configuration. This set up, although not related to critical phenomena, is 
enlightening since in this case a direct geometrical interpretation at the level of one configuration is possible, instead 
of investigating the statistical fractal properties of an ensemble. 

Firstly, we construct a finite approximation to the ID Cantor dust of 2 11 sites, with Hausdorf fractal dimension 
D f = 5/6 [im. In order to transform this fractal set into a field configuration defining a fractal measure on an 
equidistant lattice we find the minimum two-point distance of the set and using it as the lattice spacing we determine 
the size of the lattice by dividing the maximum two point distance in the set by the minimum one. In the sites of this 
new equidistant lattice that are closer to the locations of the points of the initial Cantor set, we give the field value 1, 
while in all the others we give the field value 0. Thus, we turn out with a field configuration in an equidistant lattice, 
which by construction has the property {f^°'* + \cr(x)\dx) cx C, Df in a good precision, where the reference sites xo,i 

(i = 1, 2 11 ) are obviously only those with field values 1. The averaging now is taken only on the different xo.i since 
we have only one configuration. The produced cr-field configuration is depicted in fig. [5] where the fractal property is 
clear. Note that starting from the initial 2 11 Cantor lattice, we transited in a much larger (« 2x 10 4 sites) equidistant 
one. 




FIG. 8: The a-field derived by a Cantor-like fractal, on the equidistant lattice. 

Now, we evolve this cr-configuration according to equation of motion ((3|) taking zero initial kinetic energy, and we 
focus on the evolution of (J^"' z+ '' W(x)\dx), which initially has the characteristic fractal mass dimension "0(0) = Df = 

5/6 as can be seen in the upper graph of fig. [TO] In fig. [9] we demonstrate the evolution of the mean value of the field, 
its standard deviation, and the slope ip, for potential parameters A — 1 and A = 1. We observe the same behavior 
with the case of the previous section, that is the change of the initial fractal mass dimension, and its return to this 
value each time (a) approaches zero. The explanation is the same as before without any significant new ideas. The 
only difference is that in the present case only a small subset of the entire lattice is initially described by a fractal 
measure with mass dimension 5/6. The remaining lattice sites are occupied by zero field values, which initially have 

zero contribution to (J^°' z+ '' \<j(x)\dx). However, as the system evolves, these initially zero oscillators roll towards 

the potential minimum, producing a non-fractal background, i.e a linear contribution to (J^°'"^~ \<j(x)\dx), and thus 
supplanting and suppressing the total power law behavior. On the other hand, when the mean field value returns 
to zero, the system reaches a state similar to the initial one, (zero mean field value and zero kinetic energy), the 
initially zero oscillators approach again the value a = almost in phase and reproduce a configuration with fractal 
mass dimension Df m 5/6. In fig. [10] we depict (f^"'^ \<r(x)\dx) versus £ for three times. The first corresponds to 
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FIG. 9: Mean field, standard deviation and ip (solid line) and {ip)t (dotted line) evolution for the deterministic fractal case, for 
A — 1 and A = 1, in dimensionless units. 




FIG. 10: {fx°'' + \<r{x)\dx) versus for successive times t = (initial moment), t = 2 (corresponding to the first complete 
destruction of the initial power law), and t = 4.3 ( corresponding to the first restoration of the initial fractal measure), for A = 1 
and X = 1, in dimensionless units. 



the initial moment, the second to the first complete destruction of the fractal geometry, and the third to the first 
approximate re-establishment of the initial fractal measure. 
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V. TWO COUPLED FIELDS 



It is interesting to extend the single field model of section II in the case of two coupled fields, one of them having 
a fractal profile and the other a conventional one. In this case the Lagrangian density ((T|) is extended to 

£ = \{d^a + 0„7r0"7r) - V(a, tt), (9) 

with the potential 

V(a, tt) = ~ (a 2 + tt 2 - l) 2 - Aa, (10) 

avoiding to introduce additional parameters. The fields are defined in 1 + 1 dimensions. The equations of motion 
derived from §§§ are 

a - a" + Xa 3 + A(tt 2 - l)a - A = 

tt - tt" + Att 3 + A(cr 2 - 1)tt = 0. (11) 

We are interested in studying the evolution of our system determined by the above equations, taking for the a- 
field initial conditions possessing fractal behavior as described in section II, and for the 7r-field conventional initial 
conditions, corresponding to an ideal gas at temperature To. In the next subsection we describe the generation of the 
ensemble of 1-d 7r-configurations on the lattice, using the algorithm introduced in [l4T |. 



A. Generation of thermal ^-configurations 



The unperturbed Hamiltonian for the classical scalar field theory in 1-d is 



2 

The free particle solutions for t = arc 



H =ll dx{(d t ir(x,t)) 2 + (d x ir(x,t)) 2 +mlTr(x,t) 2 }. (12) 



+ °°dk_ Akx [ + °° dk{a k + a*_ k ). 



— OO 



where LOk = \/k? + m 2 . 

Now, choosing an initial classical density distribution [3] 

p[7r,7r] = Z~\(3 ) exp{-/3 -ff[7r,7T]}, (14) 

and substitute the Hamiltonian (|12p with the free particle solutions p3p . we finally acquire 

f f +oc dk } 
p[xk,Vk\ = Z~ X (Pq) expj-A) J i^u k (xl + yl)>, (15) 

with Po = 1/2q, and where: a k = x k + iyk with Xk,yk real. In order to produce a thermal ensemble (at temperature 
To) of configurations for tt(x, 0) and tt(x, 0), we select Xk and yk from the gaussian distribution (| 1 5[) . assemble and 
then substitute in (TT5)) . The corresponding (7r(x)) profile is shown in fig.[H] All the characteristics of the 7r-ensemble 
such as the correlation function (tt(x)it(x + 8x)) — (tt(x))(tt(x + Sx)}, which turns out to be a (5-function, are consistent 
with the assumption of an ideal thermal gas. 



B. Numerical Results. Two coupled fields 



We solve equations of motion (TTT|) in 1-D 2000-site lattice following the discretization scheme ((H) , using for initial 
conditions an ensemble of 10 a and tt configurations satisfying the aforementioned requirements, that is fractal a 
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FIG. 11: The n-field on the lattice averaged over the ensemble of configurations. 
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FIG. 12: Mean field and standard deviation evolution in the two-field case for A = 1, A = 1, in dimensionless units. Mind the 
different vertical axis scale in plot c). 



and conventional 7r (corresponding to dimensionless temperature To = 1) configurations. As in the single field case 
we focus on the evolution of the (f™° + \<7(x)\dx), which defines a fractal measure, initially having the characteristic 
mass dimension 5/6, generated by the a- field. 

In figs[T2a,b we demonstrate the evolution of the mean a- field and its standard deviation, and in figures [T^b.d the 
corresponding quantities for the 7r-field, for A = 1 and A = 1. The absence of a term linear in 7r-field in the potential 
(fTTJ)) . leads to slight oscillations, due to the small initial total energy of the it- field, of the mean value (tt) around 
zero. The anharmonic character of these oscillations relies on the non- linear form of the equation of motion pip , and 
is enhanced relatively to the a case due to the large variation of the coupling term A(<7 2 — 1)tt. Note however that 
although the mean field value (it) remains small, this is not the case for its fluctuations which can become quite large. 
On the other hand, as one can see in figures [T2a,b. the corresponding a quantities do not differ significantly from the 
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single field case since 7r-field remains small. 
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FIG. 13: Time evolution of the slope 4> (solid line) and its time average {ij))t (dotted line) in the two-field case, for A — 1 and 
A = 1, in dimensionless units. 

In fig. [T3] we depict the slope ij){t) for the aforementioned evolution scenario of the a- field interacting with the 
7r-ficld. In general we observe the same phenomenon as in the single field case, that is ij)(t) becomes approximately 
equal to its initial value every time (a) approaches zero. However, the stronger non-linearity forbids (a) to become 
exactly zero but rather to obtain a small finite value, which in turn has as a consequence that ift(t) does not approach 
so closely the initial value as in the single field case. This is illustrated in figs. H2] and [T2"h . at t ~ 4 and t ~ 9, which 
correspond to the times when (<r)(t) has a local minimum. In the case of the first minimum (t ~ 4), (a) returns closer 
to zero and therefore ip(t) becomes almost 5/6, compared to the case of the second minimum at t ~ 9, when (a) 
reaches a greater value (~ 0.5) and correspondingly ip(t) deviates stronger from its initial value. However, in both 
cases the system maintains a memory of its initial fractal characteristics. 

We have evolved the system with various 7r-field initial conditions (random, constant etc) in order to check the 
generality of the results described above. It turns out that for the two-field case the re-establishment of the initial 
fractal geometry does not depend on the specific 7r-field initial conditions, unless they are fine tuned so that (a) is 
driven to large values without occasionally returning close to the region (er) w 0. 

VI. DISCUSSION AND CONCLUSIONS 

In this work we investigate the evolution of the fractal characteristics of a single or two coupled scalar fields system, 
within the framework of the general cr-model. After a relatively rapid deformation of the initial fractal geometry, we 
observe that it is being re-established almost periodically each time the mean value of the cr-field returns to zero. 
This effect is obtained using both random, as well as deterministic (Cantor-like) fractal set up. The key point for 
the occurrence of this behavior is the condition of initial equilibrium. As the system practically consists of coupled 
anharmonic oscillators, this condition is expressed through the zero initial kinetic energy for each oscillator, leading 
to a synchronous evolution of the entire system. Additionally, the fractal measure characterizing the initial state is 
based on the fact that the corresponding field configurations represent fluctuations (with a specific pattern) around 
zero. Thus, each time the oscillators pass through their turning points, associated to zero mean field value, in phase, 
the initial fractal geometry is being re-established. The above evolution, at least qualitatively, is quite general and 
robust for a very wide parametric space. 

The scenario analyzed in this paper could be extended to 3 dimensions and serve as a model in order to describe 
the evolution of the order parameter fluctuations in a critical system. In particular it could be used to explore 
the non-conventional correlations which are expected to occur during the formation of an isoscalar condensate in a 
heavy ion collision experiment. In this case one has to adapt the cr-model Lagrangian in order to describe correctly 
the characteristics of the order parameter associated with the 2nd order critical end point of the chiral QCD phase 
transition [l5| . At the phenomenological level these correlations are expressed through the fractal mass dimension 
of the cr-field configurations, determining the statistical properties of the condensate. The dynamics of the system 
is similar to the two-field case investigated above, provided that the critical system is initially at equilibrium. One 
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can study the evolution of the initial fractal characteristics of the cr-field and the possibility to leave signals at 
the detectors, supplying an indication of the phase transition. The discussion of the present work supports this 
eventuality since the time average (in order to simulate the experimental conditions) of the periodical deformation 
and re-establishment of the initial fractal geometry, will leave a signature. The specific application of the present 
work in the case of QCD phase transition has been performed in [161 ]. 
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